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I. INTRODUCTION 

Statistical moments of second order represent a key 
element of the quantum mechanical paradigm. Besides 
providing the 'language' in which the uncertainty prin- 
ciples are expressed, they serve as indicators in a num- 
ber of applications of the theory, with both applied and 
fundamental interest. In particular, second moments are 
central to the description of bosonic fields in second quan- 
tization (as is the case, for instance, in quantum optics, 
where second order coherence is characterized in terms 
of second moments) and of non-relativistic particles in 
first quantization. Such systems do, in fact, share the 
same formal description which hence extends its do- 
main to a variety of fields ranging from atomic physics to 
quantum optics, from superconductors' physics to nano- 
mechanical systems. 

During the last decade, the rise of quantum informa- 
tion science has renewed the focus on these areas because 
of their potential for coherent quantum manipulations, 
and has concomitantly brought new problems and ques- 
tions to the attention of theorists, which resulted in the 
birth of the field of "continuous variable" (CV) quantum 
information i (see Refs. 

for some literature on 
CV quantum computation, CV quantum teleportation, 
and CV quantum key distribution, respectively). A sys- 
tematic analysis of the properties of continuous variable 
quantum states inferred from the structure of their sec- 
ond moments has been thereby carried out, which lead to 
a well-established, extensive theoretical picture (see, for 
instance, [1,0j[1j[1|)- Such an analysis proved to be most 
relevant also in view of the experimental prominence of 
the class of Gaussian states [id . ITl| . which are com- 
pletely determined by their first and second moments. 
In first, seminal endeavours [1,01 j the qualitative charac- 
terization of the quantum correlations ( "entanglement" ) 
of Gaussian states of two degrees of freedom has been 
successfully achieved. Now, while very well established 
and relatively simple, this result is still often expressed 
in a non-rigourous or incomplete manner, which is prone 
to confound the unacquainted reader. Because it con- 
stitutes one of the basic building blocks on which the 



theoretical characterization of Gaussian states has been 
constructed, it seems to us extremely important for it to 
be re-derived and re-expressed in a rigourous manner and 
full detail: this is one of the motifs and central aims of 
the present paper. 

In general, the main question we will address and 
rigourously answer is the following: 

(i) What are the algebraic conditions that a 4 x 4 real 

symmetric matrix V must satisfy in order to rep- 
resent the correlation matrix of a two-mode bosonic 
system? [l^ 

Then, we shall move on to answer a closely connected 
question: 

(ii) What are the algebraic conditions to be satisfied by 

V in order to represent the correlation matrix of 
a separable (or entangled) Gaussian state of two 
bosonic modes? 

Both these questions are thoroughly answered providing 
complete sets of conditions, which are expressed in terms 
of global or local symplectic invariants. In particular, the 
set of local conditions, i.e., given in terms of local invari- 
ants, is completely new in literature. Then, by specifying 
some of these algebraic results in the case of positive- 
definite matrices, we can make a direct comparison with 
the previous work of Ref. ^ and provide a rigorous and 
correct interpretation of its seminal results. 

The paper is organized as follows. In Sec.|lT]we review 
basic notions about bosonic states, correlation matrices 
and symplectic transformations. In Sec. IIIII we present 
the mathematical tools to be used in the derivations of 
Sees. HV] and El In Sec. |TV] we provide two sets of al- 
gebraic conditions for the physical genuinity of the cor- 
relation matrix of two bosonic modes. These conditions 
are expressed in terms of global or local symplectic in- 
variants. In Sec.|V]we provide similar conditions for sep- 
arability. Next, in Section IVIi we specify some of our 
results for making a direct comparison with the previous 
achievements of Ref. 6]. Finally, Sec. IVIll is for conclu- 
sions. 
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Notice that in the paper we will denote by Ai{n,M.) 
the set of n X n real matrices. Then, we use the compact 
notation 

5(n,R) = {M e 7W(n,R) : M = M^} , (1) 

for the set of the n x n symmetric real matrices, and 

V{n, M) = {M e S{n, R) : M > 0} , (2) 

for the set of the nxn positive-definite real matrices. We 
will also consider the set (group) of proper rotations 

SO{n) = {M e 7W(n,R) : M^M = I, dctM = 1} , 

(3) 

where I is the identity matrix. 

II. BOSONIC SYSTEMS AND SYMPLECTIC 
TRANSFORMATIONS 

A. Correlation matrix of a bosonic system 

Let us consider a bosonic system composed by n 
modes, labeled by an index k. Such a system can 
be described by an infinite-dimensional Hilbert space 
H = 'S)^=i'^k and a vector of quadrature operators 
:= {qi,pi,- ■ ■ AniPn)- In particular, these operators 
satisfy the commutation relations 

[xi.Xra] = 2ir2/„ , (4) 

where Z,m = 1, - • • , 2n, and f2;„i are the entries of the 
simplectic form 

k=l ^ ^ 

An arbitrary state of the bosonic system is identified with 
a density operator p acting on the Hilbert space 7i [we 
denote by ViTi) the space of density operators acting 
on Ti]. An arbitrary density operator p e 'D{l-L) has 
an equivalent representation in a real symplectic space 
K. = /C(M^", il) called the phase space. This is a real vec- 
tor space which is spanned by the singular eigenvalues 
x^=(gi,Pi, • • • ,qn,Pn) of X"^ (representing the "contin- 
uous variables" of the system) and associated to a sym- 
plectic product u • V u-^rjv. In this space, a quantum 
state is fully described by a quasi-probability distribution 
known as the Wigner function W = W{x.). In general, 
such a function is fully characterized by the entire set of 
its statistical moments [l3|. However, in the particular 
case of Gaussian states, the Wigner function is Gaussian 
and, therefore, fully characterized by the first and second 
moments only. These two moments are also known as the 
displacement vector d := (x) and the correlation matrix 
(CM) V, whose generic entry is defined by 

Vljn ■■= ^{AxiAXm + AXm^Xl) (6) 



where Axi := xi — (xi). According to the definition of 
Eq. ([U, the CM of n bosonic modes is a real and sym- 
metric matrix in 2n dimension, i.e., V G iS(2n, R). As 
a direct consequence of Eq. Q , such a matrix must also 
satisfy the uncertainty principle 

V + . (7) 

In other words, an arbitrary V S S(2n, R) is a bona fide 
quantum CM if and only if Eq. ([7]) holds. Equivalently, 
we can introduce the set of n-mode quantum CM's to be 
defined as 

qCM{n):={VeS{2n,M.):V + in>0} . (8) 

Notice that the condition of Eq. ([7]) implies a first rele- 
vant constraint on the matrix V: 

Lemma 1 (Definite positivity of V) For every V e 
S{2n, M) satisfying V + iO > 0, one has 

V > . (9) 

Proof. Let u e R^n^ ^-^eii. < u'^Vu + iu'^Ou = 
u-^Vu because J7 is anti-symmetric. Hence V > 0. To 
prove definite positivity suppose, ad absurdum, that a 
non-trivial real vector Uq exists such that u^Vuq = 0. 
Another vector Ui € R^" such that u^Jlug ^ always 
exists as nul(r2) =0. As a consequence, one can always 
construct a set of complex vectors z = Uo + iaui, for 
a e R, such that 

< (z*)^(V + if2)z 2auf rjuo + a^uf Vui . (10) 

Values of a such that the inequality above is violated can 
always be found regardless of the values of u^JIuq and 
ufVui. This implies u-^Vu ^ for every u € R^" and, 
therefore, V > 0. ■ 

According to Lemma [TJ we then have 

qCM{n) (^V{2n,R) . (11) 

Furthermore, it is trivial to show positive-definite matri- 
ces which violate Eq. ([7]), so that we actually have 

qCM{n) (lV{2n,M) . (12) 

Notice that definite positivity is the only requirement 
for a real symmetric matrix to be a classical correlation 
matrix. 



B. Symplectic transformations 

The most general real linear transformation of the 
quadratures 

S : X — ^ x' Sx , (13) 
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must preserve Eq. (|4]) in order to be a physical operation. 
This happens when the matrix S S A4{2n,M.) preserves 
the symplectic form of Eq. ([5]), i.e., 



is a (global) symplectic invariant if 
/(V) = /(SVS^) 



(23) 



sns' =n 



(14) 



for every S £ 5p(2n,R). Then, we say that /(V) is a 
local symplectic invariant if Eq. (|23p holds for every S G 
5p(2, M) ® • • • ® Sp{2, R) tl6i] . Notice that we can extend 
the notion of symplectic invariance also to a property of 
a matrix. For instance, the definite positivity of V is a 
Sp{2n,R) := {S e 7W(2n,R) : SJ7S^ = O} , (15) global symplectic invariant since V > =^ SVS^ > 0. 



The set of all the matrices S G A^(2n,R) satisfying 
Eq. (jl4p forms the so-called real symplectic group 



whose elements are called symplectic or canonical trans- 
formations. As a consequence, the most general real lin- 
ear transformation in phase space S : x — > x' := Sx 
must be symplectic. Its action on the Wigner function 
is simply given by W{k) — > W{S~^x), so that the dis- 
placement is linearly modified while the CM is trans- 
formed according to the congruence 



V 



svs' 



(16) 



Symplectic transformations are very important since ev- 
ery S acting in the phase space /C corresponds to a Gaus- 
sian unitary U{S) acting on the Hilbert space Ti, i.e., a 
unitary operator preserving the Gaussian statistics of the 
quantum states. These unitaries are the ones generated 
by bilinear Hamiltonians and can always be decomposed 
into sing le-mode squeezers and multi-mode interferome- 
ters flQlllsll. In particular, local symplectic transforma- 
tions 



S=0Sfce5p(2,] 



>Sp(2,] 



k=l 



correspond to local Gaussian unitaries 



U{S) = (g) Uk 

k=l 



(17) 



(18) 



Local symplectic transformations can always be decom- 
posed as products of local rotations and local squeezings. 
In fact, thanks to the following characterization 



5p(2, 



{S e 7W(2,R) : detS = 1} , (19) 



we have that every S G 5p(2,M) can be expressed as a 
product of proper rotations 



and squeezing matrices 
S(0 := 



sm if 

C0Sly9 



- cos if 
simp 



^1/2 







^-1/2 



(20) 



(21) 



Besides, it is also important to identify which quan- 
tities of a CM are preserved under the application of 
symplectic transformations. In general, for a given CM 
V, we say that a functional 



/ : V -> /(V) G 



(22) 



III. SYMPLECTIC ANALYSIS 

Here, we review some basic tools that can be used for 
the symplectic manipulation of the CM's. In particular, 
the central tool in this trade is Williamson's theorem [l3] , 
which ensures the possibility of carrying out the symplec- 
tic diagonalization of real matrices in even dimension un- 
der the definite positivity constraint (as in the case of the 
CM's). 

Lemma 2 (Williamson's theorem) For every V G 
7^(2n,R), there exists a symplectic matrix S G iSp(2n,R) 
such that 



SVS' 



/-I 



v 



1^1 



\ 



W > 



(24) 



where the n positive quantities {i>i, • • • , i>„} are called the 
"symplectic eigenvalues" ofV, and the diagonal matrix 
W is called the "Williamson form" (or "normal form") 
ofV. 

The symplectic spectrum {vi,--- ,Vn} can be com- 
puted as the standard eigenspectrum of the matrix |ir2V| 
where the modulus must be understood in the operato- 
rial sense [l^ . The corresponding Williamson form W is 
unique up to a permutation of the symplectic spectrum 
(i.e., of the bosonic modes). Fixing this permutation, the 
diagonalizing symplectic matrix S of Eq. is defined 
up to local rotations ©JJ^iRfc with Rfe G 50(2). For the 
sake of completeness, we report in Appendix [Xl a simple 
proof of Williamson's theorem, originally presented in 
Ref. [l^ (see also Ref. [HI). Using this proof, we show 
in Appendix [B] an algorithm which finds the diagonaliz- 
ing symplectic matrix S of Eq. ((24| . This algorithm is 
not the fastest but it can be helpful in studying prob- 
lems like the optimal discrimination of Gaussian states 
20] and the Quantum Illumination [21]. 

Let us consider the case of a 4 x 4 positive-definite 
matrix V G 7^(4, R), as in the case of CM's describing 
two bosonic modes. This matrix can be expressed in the 
blockform 



V = 



A C 
B 



(25) 
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where A,B e 5(2, K) and C G A^(2,M). In this case the 
symplectic spectrum {1^1,1^2} can be com- 

puted via the simple formula [13] 

where 

A(V) := detA + detB + 2dctC . (27) 

Here, the quantities dot A, dot B and det C are local sym- 
plectic invariants, while det V and A(V) are global sym- 
plectic invariants, which can also be written as 

det V ^iyW+ , A(V) = + i^l . (28) 

Another important tool in the symplectic analysis is 
the reduction to standard form by local symplectic trans- 
formations @, 0| • In general, such a reduction holds for 
symmetric matrices V G 5(4, M) with positive diagonal 
blocks, as we easily show in the following. In particular, it 
can be applied to positive-definite matrices V G 7^(4, R) 
and, therefore, to CMs V G qCMi2). 



where 

A' : = RieA) (Sa a S^) R{9Af , (35) 

B' : = II{9b) (Sb B S^) R(^b)^ , (36) 

C : = R(0a) (Sa C S^) Riesf . (37) 

Since A,B G 7^(2, M), we can apply Williamson's theo- 
rem. This means that we can choose Sa and such 
that 

A' = R(6Ia) al R.{eAf = al , (38) 
B' = R(6is) 61 R{0Bf 61 , (39) 

where a (b) is the symplectic eigenvalue of A (B) while 
the angle 6a (Ob) is arbitrary. Since the pair {6a, Ob} 
can be chosen freely, we can always choose a pair {0a, 0b} 
in Eq. ([37|) such that C = diag(c+,c_) [2^. As a con- 
sequence, Eq. is globally equal to Eq. ([50]) . Finally, 
since the transformation of Eq. ([55]) is local and symplec- 
tic, all the determinants relative to the blocks and the 
global matrix are preserved, so that Eqs. (PT|) and ([5^ 
are trivially implied. ■ 



Lemma 3 (Standard Form) For every 

V= (^^T b) e5(4,R) , (29) 

with A, B > 0, there exists some S G 5p(2, M) 5p(2, M) 
such that 




where the real parameters a, b, €+, c_ satisfy 

det A = , det B = 6^ , det C c+c_ , (31) 

and 

detV = dctV^ = {ab - cl){ab - c^_) . (32) 



IV. GENUINENESS OF A TWO-MODE 
CORRELATION MATRIX 

By applying the symplectic tools of the previous 
Sec. IIII[ we can now derive very simple algebraic con- 
ditions for characterizing the genuineness of a two-mode 
CM. In other words, starting from a generic 4x4 real 
and symmetric matrix [V G 5(4, M)], we give the alge- 
braic conditions that such a matrix must satisfy in order 
to represent the CM of two bosonic modes, i.e., a bona 
fide two-mode quantum CM [V G qCM{2)]. As a con- 
sequence of Williamson's theorem, we have the following 
algebraic conditions in terms of global symplectic invari- 
ants 0. 

Theorem 4 An arbitrary V G 5(4, M) is a quantum CM 
if and only if it satisfies 

V > , iy->l , (40) 



Proof. Let us consider a pair of single-mode symplectic 
transformations 8^1,83 G 5p(2,R) and a pair of single- 
mode proper rotations II{0a),R{9b) G 50(2). By ap- 
plying the local symplectic transformation 

8 = R(^^^)8A©R(0s)8s (33) 
to the matrix V, we get 

8V8^ = ( (^T ) , (34) 



or, equivalently, 

V > , det V > 1 , A(V) < 1 + det V . (41) 

Proof. For every V G V{A,K), the application of 
Williamson's theorem to Eq. ([7]) gives V -I- ifi > ■^f=^ 
V- > 1 (recalling that is the smallest symplectic 
eigenvalue). Since V -I- zil > => V > (see Lemma 
1), we can write V iJl > (V > A > I) 
which proves the bona fide condition of Eq. for a 
generic V G 5(4, M). Under the definite positivity as- 
sumption V > 0, one can also use Eq. to prove the 
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equivalences 

V- > 1 
max 



=^ A(V) - 2 > ^A(V)2 -4dctV 
;{2, 2\/det V} < A(V) < 1 + det V 
det V > 1 

2VdetV< A(V) < 1 + dctV . ^ ' 



According to Eq. ^ the condition 2VdetV < A(V) 
in Eq. ((42|) corresponds to 2i'^i'^ < ^ which is 
trivially satisfied. Then, we have 



dot V > 1 
A(V) < 1 + detV 



(43) 



which states the equivalence between to Eqs. (j40|) 
and (|4ip . where the underlying assumption V > is 
also shown. ■ 

Notice that, crucially, Williamson's theorem could be 
applied to V because of its definite positivity, which thus 
implies the existence of well-defined symplectic eigenval- 
ues [2^. The condition i^- > 1 alone is therefore not, 
by itself, fully equivalent to the uncertainty principle 
V + in > 0, unless definite positivity is also assumed. 
The essential role of the prescription V > 0, often ne- 
glected in the literature, is especially clear in the formu- 
lation of Eq. (|¥T|l : in fact, the other two inequalities in 
Eq. (|4ip only depend on the squared symplectic eigen- 
values and cannot thus distinguish between positive and 
negative eigenvalues. 

Besides the algebraic requirements of the previous the- 
orem, we can derive an alternative set of conditions by 
applying the reduction to standard form. These new con- 
ditions are expressed in terms of local symplectic invari- 
ants. 



Theorem 5 An arbitrary 
A 



e 5(4,] 



(44) 



is a quantum CM if and only if it satisfies 

A,B > , 
A(V) < 1 + detV , 



(45) 
(46) 

2%/det AdctB + (det C)^ < det V + det A det B . (47) 

Proof. Under the assumption A, B > 0, the matrix V 
can be reduced to the standard form of Eq. ([30| via 
a local symplectic transformation S. Since 
the Heisenberg principle can be written in the equivalent 
form fH 



V + in>0 



in> 



(48) 



Since the matrix + id is Hcrmitian, its four eigen- 
values A]{], At, A^, AI arc real. It is then easy to show 



that 



2A;^ = a + b+ Jfi± 2^ 



2A^ = a + 6 - -y//7±2V^ 



(49) 
(50) 



where 



and 



4+(a-5)2 + 2(c^ + ct)>4, (51) 



V ■= 4(a - hf + (c+ + c_)2[4 + (c+ - c_)2] > . (52) 
Since is the minimum eigenvalue, we have that 

+ irj > AT > ^ 



(a + 6)2 > ^ + 2V^ 
a-f 6> 



/^ + 2Vi^> 

[(a + 5)2 - > 41/ 
(a + 6)2 - ^ > 
a + 6 > . 

(53) 



Last condition a + 6 > in Eq. ([55)1 is trivially included 
in A > and B > which gives a > and 6 > (for 
congruence with the diagonal matrices al and 51). The 
other two conditions 



and 



[(a + 6)2 - ^lf > 



(54) 



(55) 



can be recast in terms of the local symplectic invariants. 
In fact, by inserting Eqs. ([5T|) and ((5^ in Eq. ([5i|) . we 
get 



0^ + 6^ + 2c+c_ < [ah - c\){ab - c^) + 1 



(56) 



which is equivalent to Eq. (HS)) by using Eq. ([5^ and 
A(V) = A(V^) = a2 + 62 + 2c+c_. Finally, by inserting 
Eq. dnil) in Eq. ([551), we get 2a5 - - c2_ > 2 which is 
equivalent to 



2a''6 



2>2 



a6(c^ + c^) > 2a5 , 



(57) 



since a6 > 0. In terms of local symplectic invariants, last 
inequality is equal to 



2 det A det B - /4 > 2\/det A det B , (58) 



where 



/4 := Tr(Aa;Ca;Ba;C'^a;) = a6(4 + ct) (59) 

is another local symplectic invariant. In fact, the quan- 
tity /4 is connected to the other local symplectic invari- 
ants by 



det V = det A det B + (det C) - 1^ 



(60) 



which holds for every V g 5(4, M). Using Eq. (pD|) in 
Eq. diSD, we then get Eq. (gT]). ■ 
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V. SEPARABILITY OF A TWO-MODE 
CORRELATION MATRIX 

Few years ago, Ref. ^] showed how to extend the par- 
tial transposition and the Peres entanglement criterion 
[2^ to bipartite bosonic systems. In fact, partial trans- 
position PT : pAB — > PAB corresponds in phase space 
to a "local time reversal" which inverts the momentum 
of only one of two subsystems. This means that we have 
the following transformation for the Wigner function 



PT : VF(x) — > T^(x) VF(Ax) , 



where 



A := 



1 



-1 



(61) 



(62) 



For the corresponding CM V e gCA1(2), the PT trans- 
formation is given by 



PT : V 



V := AVA 



(63) 



where the partially transposed matrix V belongs to 
7^(4, K) [Hi but not necessarily to qCM{2). By writ- 
ing V in the blockform of Eq. (pS)) . one easily checks 
that the action of the PT transformation A reduces to 
the following sign flip 



det C ^ - det C , 



(64) 



at the level of the local symplectic invariants. As a con- 
sequence, the positive-definite matrix V has 



A(V) = detA-f detB-2detC 
and symplectic eigenvalues 



A(V) , (65) 



A(V)±i/A(V)2~4detV 



(66) 



Once the PT transformation has been extended, also the 
Peres criterion can be consequently extended via the log- 
ical implication 



PAB separable =^ pab G I?(H) 
=^ V e qCM{2) , 



(67) 



which becomes an equivalence for Gaussian states under 
1 X n mode bipartitions 0, Q . 

Theorem 6 (Separability) Let us consider a Gaus- 
sian state PAB with CM V G qCM(2). Then, pab is 
separable if and only if 



V e qCM[2) , 



or, equivalently, 
or, equivalently, 



A(V) < 1-hdetV 



(68) 
(69) 
(70) 



Proof. The proof of Eq. l|68|) follows exactly the same 
steps of the one in Ref. |6|, where the P-representation 
is exploited (see Ref. [23 for recent connections be- 
tween P-representation and separability.) In order to 
prove Eqs. ([69]) and ([70]), let us apply Theorem [4] to the 
positive-definite matrix V G 7^(4, M). Then, we get 



V e qCM{2) 



> 1 



detV >1 
A(V) < 1 + detV 

(71) 

where Eq. (I69p is trivially proven. Now, since 
detV = dot (AVA) = detV > 1, the first condition in 
Eq. ([7T|) is always satisfied and, therefore, the separabil- 
ity condition is reduced to Eq. ([70)1 . ■ 

Let us now derive the algebraic conditions that a 
generic symmetric matrix must satisfy to represent the 
CM of a separable or entangled Gaussian state. The 
following corollary gives an easy recipe to check if a sym- 
metric matrix is a good or bad candidate for this aim. 



Corollary 7 An arbitrary 
V = 



A C 
B 



e 5(4, M) 



(72) 



represents the CM of a separable Gaussian state if and 
only if it satisfies 



V > , iy->l, > 1 , 
or, equivalently, 

V > , dctV > 1 , r(V) < 1 + dctV 
or, equivalently, 

A,B >0 , r(V) < 1 + detV , 



(73) 



(74) 



(75) 



2%/det AdetB + (det C)^ < det V det A det B , (76) 

where T(V) := det A -I- det B -I- 2 |det C|. Instead, it repre- 
sents the CM of an entangled Gaussian state if and only 
if it satisfies 



V > , iy->l, 
or, equivalently, 

V > , detV > 1 , A(V) < 1 + detV < A(V) 
or, equivalently, 

A,B > , A(V) < 1 + detV < A(V) , 



2 Vdet A det B (det C)^ < det V det A det B 



(77) 

(78) 

(79) 
(80) 



Proof. In order to represent the CM of a separable 
Gaussian state, the symmetric matrix V e 5(4, R) must 
simultaneously satisfy 



V e qCM{2) , V e qCM{2) 



(81) 
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The bona fide condition V G qCM (2) is equivalently ex- 
pressed by the conditions of Eqs. and (jH)) in The- 
orem HI Then, for every V G gCA4(2), the separability 
condition V G (7C7W(2) is equivalent to Eqs. and (jTO]) 
in Theorem [H By combining Eq. (|40|) with Eq. ([69]) . 
and Eq. (gl]) with Eq. jTOl), one easily gets Eqs. jTSl) 
and dill), where max{A(V), A(V)} < 1 detV 
r(V) < 1 -I- det V. According to Theorem [6l for every 
V G qCM.{2) the entanglement condition is expressed by 
< 1 or, equivalently, by A(V) > 1 -I- det V. Then, it 
is trivial to derive the corresponding Eqs. (|77p and ([75)1 . 
The proof of Eqs. ((75ll76l) and ((79ll80)) is the same as be- 
fore except that now we have to combine the Eqs. (j45|) . 

and (gT]) of Theorem E] with Eq. for the separa- 
bility and with A(V) > 1+det V for the entanglement. ■ 



In terms of local symplectic invariants, last relation can 
be written as in Eq. ([55)) . where Eq. ([5^]) has been also 
used. ■ 

Notice that Eq. ([55]) corresponds to the Eq. (17) of 
Ref. ^Qi] , up to notation factors 28] . In Ref. [6,] , this condi- 
tion is incorrectly claimed to be equivalent to the Heisen- 
berg principle V + i^l > (this equivalence is claimed 
under the positivity contraint V > 0, which is a suffi- 
cient condition for the reduction to standard form used 
in the proof of Ref. Q)- In order to have a full equiv- 
alence with the Heisenberg principle V -f ifi > 0, the 
supplementary condition of Eq. (j84p is mandatory. It is 
indeed rather simple to construct a positive-definite ma- 
trix V G P(4,R) which satisfies Eq. jHS]) but violates 
V-l-iJl > 0. As an example, let us consider the following 
real and symmetric matrix 



VI. RELATION WITH THE PREVIOUS WORK 
BY SIMON 

In order to make a direct comparison with the previous 
work by Simon Q , we have to specify some of our results, 
given for arbitrary symmetric matrices V G 5(4, K), to 
the case of positive-definite matrices, i.e., V G P(4,R). 
As an immediate consequence of Theorem |4l we have the 
following result. 

Corollary 8 An arbitrary V G 7^(4, M) is a two-mode 
quantum CMY G qCM{2), i.e., V + iJl > 0, if and only 

i/_ > 1 , (82) 

or, equivalently, 

detV > 1 , A(V) < 1 + detV , (83) 
or, equivalently, 

detV>l, (84) 

det A det B (1 - det C)^ - /4 > dot A + dot B , (85) 

where h Tr(Aa;Cu;Bu;C'^u;). 

Proof. By applying Theorem [4] under the assumption 
V > 0, one trivially derives the equivalent conditions in 
Eqs. ^ and In order to prove Eqs. ^ and ([55)l . 
let us reduce the positive-definite matrix V to its stan- 
dard form of Eq. ([50)1 . Under local symplectic transfor- 
mations, we then have the equivalence 

A(V) < 1 + dctV^ 
a2 + 6^ + 2c+c_ < 1 + {ab - cl){ab - d) . (86) 

Note that Eq. ([86]) can be equivalently written as 



/ l + 4x -l + Ax \ 
, - 1 I 1 -I- 4a; -ix 1 
^ 2 -1 + 4a; 1 + 4a; 

V -4a; 1 -I- 4a; / 

(88) 

which is positive-definite for every a; > 0. It is easy to 
verify that the Hermitian matrix V+ifJ has the following 
real eigenvalues 

A± = i(l + 8a;± Vl7-16a; + 64a;2) , (89) 

d± = i(3 + 8a;± Vl7- 16a; + 64a;2) . (90) 

Since A_ is the minimal eigenvalue, the Heisenberg prin- 
ciple V -|- ill > is equivalent to A_ > 0, which gives 

a; > 1/2 . (91) 

Now, let us explicitly compute Eqs. ([84| and ([85| . It is 
easy to show that Eq. ([84]) is equivalent to 

x(8a; + 1) > 1 ^ x> (%/33 - 1)/16 ~ 0.3 , (92) 

while Eq. (I85p (Simon's genuineness condition) is equiv- 
alent to 

1 -I- a;(8a; -5)>0 4=» 0<a;<ioRa;>i. (93) 

2 8 2 

From Eq. ([93]) . one can see that Simon's condition 
alone does not exclude the matrices V(a;) for < a; < 
1/8, which are clearly unphysical since they violate the 
Heisenberg condition of Eq. (pij) . A complete equiva- 
lence with Eq. (PTjl is retrieved by coupling Eq. ([M)) with 
Eq. ([m, where the latter equation excludes the non- 
physical region < a; < 1/8. 

This imprecision in Simon's work leads to a common 
misunderstanding of the subsequent separability condi- 
tion [Eq. (19) of Ref 0], which in our notation corre- 
sponds to 



a^b^ -f (1 - c+c-f - ab{cl + c^_)>a^ + b^ 



(87) 



detAdetB-t-(l - |detC|)^-/4 >detA + detB . (94) 
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In this condition, the Heisenberg principle is erroneously 
claimed to be included (in fact, it is only partially 
included). Hence, Simon's separability condition of 
Eq. ([M]) is actually vahd only if V € qCM{2), i.e., the 
positive-definite matrix V e 7^(4, R) is already known to 
be a bona fide quantum CM. In other words, the sepa- 
rability criterion of Eq. ([M]) must be tested on positive- 
definite matrices which are already known to describe 
the second statistical moments of a physical quantum 
state. However, under this assumption of physicality, 
Simon's separability criterion of Eq. (|94p displays a re- 
dundant modulus and must be simplified to 

det A det B + (1 + det C)^ - /4 > dot A + det B . (95) 

Criterion 9 (Separability) Let us consider a two- 
mode quantum state pab having quantum CM 

V= (^^T eqCM{2) . (96) 

The separability of pab implies 

A(V)<l + detV, (97) 

or, equivalently, 

det A det B + (1 + det C)^ - /4 > det A + det B . (98) 

In particular, if pab is a Gaussian state, then it is sep- 
arable if and only if Eq. |g7| ) [or Eq. i98\) _l holds. 

Proof. The proof is straightforward. Since V is a quan- 
tum CM, it is positive-definite and satisfies the condition 
detV > 1. Now, suppose that the corresponding two- 
mode state Pab is separable. Then we have 

PAB separable =4> pab &T>{H) 

V e qCM{2) + in>0 . (99) 

By applying Corollary [5] to the positive-definite matrix 
V = AVA, we have 

V + iJ7 > ^ detV > 1 , A(V) < 1 + detV . (100) 

Since det V — det V, we have that det V > 1 is auto- 
matically satisfied in the previous Eq. (jlOOp . Then, we 
get 

PAB separable =^ A(V) < 1 det V , (101) 

where A(V) := A(V) is defined in Eq. Us- 
ing Eq. (|60p . one easily proves the equivalence between 
Eqs. ([57)1 and ([55)1 . Finally, the full equivalence which 
holds for Gaussian pab is a direct application of Theo- 
rem |6l H 

This criterion represents a simplification of Simon's 
separability criterion. Now, it is important to notice that 
the separability criterion becomes a bit more involved 



when arbitrary positive-definite matrices V G 7^(4, R) 
are considered, without any other a priori assumption. 
For a generic V G 7-'(4, R), both the Heisenberg principle 
(V-t-zil > 0) and the separability property (V-t-zfJ > 0) 
must be explicitly considered and combined together, 
in order to get a complete set of algebraic conditions. 
Thanks to these conditions, one easily checks when a 
positive-definite matrix V € 7^(4, M) can represent the 
quantum CM of a Gaussian state pab which is separa- 
ble or entangled. By applying Corollary [71 we get the 
following criterion for positive-definite matrices. 

Criterion 10 An arbitrary V e 7^(4, R) represents the 
CM of a separable Caussian state if and only if 

detV>l, (102) 

r(V)<l + dctV, (103) 

or, equivalently, 

detV > 1 , (104) 

det A det B (1- |detC|)^ ~ h> det A + det B . 

(105) 

Instead, it represents the CM of an entangled Caussian 
state if and only if 

detV > 1 , (106) 
A(V) < 1 + detV < A(V) , (107) 

or, equivalently, 

detV>l, (108) 
(1 + det C)^ < det A + det B - det A det B -\- h 

<(l-detC)^ (109) 

The proof is a trivial application of Corollary [7l to- 
gether with Eq. ([50]) . used to state the equivalences be- 
tween Eq. pug)) and Eq. pll5)) . and between Eq. pIT7|) 
and Eq. (|109p . According to Eq. (|109p . positive-definite 
matrices with det C > can only be associated to separa- 
ble Gaussian states @]. Notice that the original Simon's 
separability criterion, i.e., Eq. (jl05p . must be coupled 
with the mandatory condition of Eq. (|104p in order to in- 
vestigate correctly the separability properties of a generic 
positive-definite matrix. 

VII. SUMMARY 

In Theorem |4l we have re-derived and explicitly stated 
all the precise algebraic conditions a symmetric matrix 
must satisfy to represent the CM of a two-mode bosonic 
(or canonical) quantum system, including the (critical 
and often neglected) definite positivity condition. Such 
conditions are expressed in terms of global symplectic in- 
variants. In Theorem[51 we have derived a new and alter- 
native set of conditions, which are expressed in terms of 
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local symplectic invariants. In these local conditions the 
positivity check is restricted to the submatrices A and B. 
Finally, in Theorem [6l the necessary and sufficient condi- 
tion for the separability of two-mode Gaussian states has 
been reviewed and cast in a compact form. We should 
stress that such a condition is valid only under the as- 
sumption that physicality is also met [V e qC A4{2)]. In 
Corollary [71 both the physicality and separability have 
been explicitly taken into account. Then, we have de- 
rived a complete set of (global or local) conditions that a 
generic symmetric matrix must satisfy in order to repre- 
sent the CM of a separable (or entangled) Gaussian state 
of two bosonic modes. In Section IVll some of our results 
have been specified for positive-definite matrices and a 
comparison with the previous results by Simon has been 
thoroughly presented. 

The rigourous agreement with all the conditions here 
considered should constitute a constant reference in both 
the theoretical practice and the analysis of experimental 
data involving quantum systems of two canonical degrees 
of freedom. 



are antisymmetric (because V and W are symmetric, 
while n is antisymmetric). In particular, we have 

k=l ^ *■ 

Now the existence of R in Eq. (|A3[) is assured by the fol- 
lowing theorem on the block-diagonalization of real anti- 
symmetric matrices (specialized to even dimensions) (29| 

Theorem 11 For every A = — A^ e M{2n,R), there 
exists a (unique) O G SO{2n) such that 

n 

OAO^ = 0afea;:=A , (A6) 

fc=i 

where the (unique) block diagonal form A has Ok > 0. 

APPENDIX B: FINDING THE DIAGONALIZING 
SYMPLECTIC MATRIX 
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APPENDIX A: SIMPLE PROOF OF 
WILLIAMSON'S THEOREM 

Let us construct the diagonalizing symplectic accord- 
ing to the decomposition 

S = W^/^RV-i/^ ^ (A^) 
with a suitable R G SO{2n). In fact 

SVS^ = (W1/2RV-1/2) V (V-i/'R^Wi/2) 
= W^/^j^ I R'^W^/^ ^ ^1/2 J ^1/2 ^ ^ ^ ^^2) 

Notice that Eq. (|Aip is well-defined since V and W are 
positive-definite (therefore, non-singular). However, the 
rotation R in Eq. (|Aip is not arbitrary but must be cho- 
sen in order to make S symplectic. 

Let us apply Eq. (jAl[) to the symplectic condition 
SflS^ = n. Then, we have 

(W1/2RV-1/2) n (v-i/'R^wi/2) = 

^ R (V-^/'OV-i/2) R^ = W-i/2f2w-i/2^ 

^ RXR^ = Y , (A3) 

where 

X := V"i/2f2V-i/2 ^ Y W^^/^^jw^^/^ (A4) 



Let us show a possible procedure for deriving the 
proper rotation O that block-diagonalizes a generic an- 
tisymmetric matrix A as in Theorem 1111 We can eas- 
ily prove the following connection between the block- 
diagonalization of A and its unitary diagonalization 

Theorem 12 The proper rotation O performing the 
block-diagonalization of Eq. HA0) is given by 

o = ru^ , (Bi) 

where 

and U is an arbitrary unitary performing the diagonal- 
ization of A, i.e., 

UtAU = 0zafe ( ^ j := Ai, . (B3) 
fe=i ^ ^ 

Proof. First, let us prove how A can be transformed 
into the diagonal form Ajj of Eq. (|B3|) by a unitary ma- 
trix. From Eq. (jB2[) . we have that 

7'f (flfcw) 7 = iflfc ^""^ -^^ . (B4) 
As a consequence, by applying F to Eq. (|A6[) . we get 
r^OAO^r = Ft AT = ^iaJ'^ 1 ) = ' (B^) 

In other words, there exists a unitary O-'^F that diago- 
nalizes A according to Eq. (|B3[) . 



10 



Then, let us prove that, for every unitary U diagonal- 
izing A according to Eq. (|B3p . we can write Eq. (|B1[) 
where O performs the block-diagonalization of Eq. (jA6|l . 
For proving this, let us consider the orthonormal eigen- 
vectors {Ui, • • • , U2„} of A 



Aui = — miUi , Au2 = iaiU2 , 



(B6) 



Au2„-i = -ia„U2„-i , Au2„ = +m„U2„ , (B7) 

more compactly denoted by {u2k-i,Vi2k}k=i with 

Au2fc_i = -iafcU2fc-i , Au2fc = iakU2k ■ (B8) 

These vectors are unique up to phase factors <p := 
{(pi, • • • , </32n}, i-e., up to the replacements 



U2fc-1 U2k-ie"f""'-^ , U2k 



U2fce 



(B9) 



This means that, for every choice of we have an equiv- 
alent unitary matrix U — U(<p) in the diagonalization of 
A. Now, by conjugating Eq. (|B8p . one easily checks that 



U2fc-1 



^2k ■ 



(BIO) 



As a consequence, the most general unitary matrix that 
diagonalizes A has the specific form 



U = ( U2 • • • U^„ U2„ 



(Bll) 



Let us explicitly compute the matrix product FU'''. By 
applying 



/ i —i 
1 1 



V 







I —I 
1 1 J 



(B12) 



to the conjugate matrix 



ut 



^2n 
\4nJ 



U2A U2.2 



U2,2n~l U2,2n \ 



U 



2,1 



U 



2,2 



^2,2n-l 



*2,2n 



U2n,l U2n,2 
\ "2n,l "2n,2 

one explicitly gets 



'W2n.2ri-l U2n,2r. 



2n,2n-l 



(B13) 



!n,2n / 



ru+ 



V2 



( ai,i ai,2 

/9l,2 



an, 2 



OL\.2n-\ ai,2n \ 
Pl,2n-1 Pl,2n 



Otn,2n-l Oln,2n 
Pn,2n-1 Pn,2n I 



(B14) 



where 



ttfej := -2Im(M2fe,i) , /3fe,j := 2Re(u2fe,i) 



(B15) 



From Eqs. (|B14|1 and (|B15p . we have that TU''' is real 
for every choice of U in Eq. (|B3p . i.e., for every choice 
of the phases in the corresponding eigenvectors. More 
strongly, we have FU''' e 50(2^1) (since FU''' real implies 
FlJt orthogonal with det = +1). Then, from Eq. (jB3|l . 
we easily get 



FU+AUF^ = FAnF^ = A 



(B16) 



In conclusion, for every diagonalizing unitary U, the 
proper rotation FU^ corresponds to the unique proper 
rotation O that performs the block-diagonalization of 

Eq. (Hg. ■ 

Both Theorem [11] and Theorem [12] can be applied to 
the Eq. ^M^, by setting O = R, A = X and A = Y. 
These theorems allow to reduce the computation of the 
rotation R in Eq. (|A3|1 to a unitary diagonalization. In 
fact, we have just to find a unitary U that diagonalizes 
X, i.e.. 



k=l 



and then construct 



R = FU+ 



(B17) 



(B18) 



Once that we have R, we use Eq. (|Aip to get the sym- 
plectic S. Here is the complete algorithm: 

1. Find the symplectic spectrum of V, i.e., its 
Williamson form W 

2. Compute the matrices W^/^ (immediate) and 

(needs orthogonal diagonalization) 

3. Construct the matrix X ;= V^^/^fJV"^/^ 

4. Find the eigenvectors of X and construct the cor- 
responding unitary U 

5. Compute R = FUt 

6. Compute S = W^/^py-i/a 

By construction, this algorithm reduces the determina- 
tion of S to unitary diagonalizations. Actually, this task 
can be achieved via faster methods when the symplectic 
spectrum is non-degenerate. In general, the determina- 
tion of S is equivalent to the construction of a symplectic 
basis [30] . 
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